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Qh' Abstract 
D ■ 

I For = 1 supersymmetric quantum electrodynamics, regularized by higher 

derivatives, a method for summation of all Feynman diagrams defining the /3- 
function is presented. Using this method we prove that the /3-function is given 
by an integral of a total derivative, which can be easily calculated. It is shown that 
' surviving terms give the exact NSVZ /3-function. The results are compared with 

. the explicit three-loop calculation. 
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2 '■ 1 Introduction. 



Quantum corrections in supersymmetric theories (and, in particular, the /3-function) 

exact /3-function for = 1 sn 

3C2-T{R) + C{R),hi{(^)/r' 



^ , were studied for a long time. The exact /3-function for = 1 supersymmetric theories 



/3(a) 



27r(l - C2a/27r) 

was found in pj. Certainly, a /3-function in supersymmetric theories was also calculated 
explicitly in the lowest loops. Most calculations were made with the dimensional reduc- 
tion |2j in MS'-scheme [3]. The NSVZ /3-function agrees with these calculations in the 
one- and two-loop approximations. In order to obtain NSVZ /3-function in higher loops 
it is necessary to perform a special redefinition of the coupling constant The possi- 
bility of such a redefinition is very nontrivial [5J. However, it is well known [6] that the 
dimensional reduction is not self-consistent. (The dimensional regularization [7] breaks 
the supersymmetry and is not convenient for calculations in supersymmetric theories.) 
Ways allowing to avoid such problems are discussed in the literature [8]. 

Other regularizations are also sometimes applied for calculations in supersymmetric 
theories. For example, in Ref. [9] a two-loop /3-function of the N = 1 supersymmet- 
ric Yang-Mills theory was calculated with the differential renormalization ^U\. Some 
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calculations were made with the higher covariant derivative regularization, proposed in 
[TT] . which was generalized to the supersymmetric case in Ref. [12] (another variant 
was proposed in Ref. [I3]). Usually integrals arising with the higher covariant deriva- 
tive regularization can not be calculated analytically in higher loops. That is why this 
regularization was applied for explicit calculations rather rarely. In particular, the first 
calculation of quantum corrections for the (non-supersymmetric) Yang-Mills theory was 
made in Ref. [H]. Taking into account corrections, made in subsequent papers [15], the 
result for the /3-function appeared to be the same as the well-known result obtained with 
the dimensional regularization [16]. In principle, it is possible to prove that in the one- 
loop approximation calculations with the higher covariant derivative regularization always 
agree with the results of calculations with the dimensional regularization [ITj. Some cal- 
culations in the one- and two-loop approximations were made for various theories [TK| [T^ 
with a variant of the higher covariant derivative regularization, proposed in [20]. The 
structure of the corresponding integrals was discussed in Ref. [T9] . 

The three-loop /3-function for the = 1 SQED, regularized by higher derivatives, 
was calculated in ^21j. This calculation shows that the integrals defining the /3-function 
are integrals of total derivatives. A two-loop calculation made with the dimensional 
reduction in the N = 1 SQED and revealing a similar feature was presented in Ref. 
[22], where the factorization of integrands into total derivatives is explained in all loops 
using a special technique, based on the covariant Feynman rules in the background field 
method [23]. This factorization allows to calculate one of the loop integrals analytically. 
As a consequence, (in the A^ = 1 SQED) integrals defining the /3-function are reduced to 
integrals defining the anomalous dimension, producing NSVZ /3-function. In particular, 
with the higher derivative regularization it is not necessary to perform a redefinition of 
the coupling constant. The factorization of integrands into total derivatives seems to 
be a general feature of supersymmetric theories. In the two-loop approximation it was 
verified for a general renormalizable A^ = 1 supersymmetric theory, regularized by higher 
covariant derivatives, in |24] . 

Using a usual supergraph technique [251 126] an attempt to prove the factorization of 
integrands into total derivatives was made in Ref. [27], where a solution of the Ward 
identity was substituted into the Schwinger-Dyson equations. This allows to present a 
/3-function as a sum of two contributions: The first one (calculated in Ref. [27] exactly 
to all orders) is given by an integral of a total derivative and is expressed in terms of the 
two-point function of the matter superfields. This contribution gives the exact NSVZ /5- 
function. The second contribution is essential starting from the three-loop approximation 
and can be expressed in terms of a transversal part of a certain Green function. Explicit 
three-loop calculations [21] show that the second contribution to the /3-function is also 
given by an integral of a total derivative and is equal to 0. It was conjectured [27] that 
this takes place in all loops. A partial four-loop verification of this statement were made 
in Ref. [28]. The factorization of the additional contribution into an integral of a total 
derivative was qualitatively explained in [29], using a technique proposed in Ref. [30] . In 
this paper we formulate these arguments in a rigorous form and prove that all integrals in 
the A^ = 1 SQED defining the /3-function are integrals of total derivatives. Taking these 
integrals we obtain the exact NSVZ /3-function without a redefinition of the coupling 
constant. The results are compared with the explicit three-loop calculation, made in [2T] . 

The paper is organized as follows: 

In Sec. [2] we recall basic information about the A^ = 1 supersymmetric electrodynamics 
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and its regularization by higher derivatives. In Sec. [3] we perform the integration over 
the matter superfields in the generating functional and introduce some notation. Then it 
is possible to construct formal expressions for a /3-function and an anomalous dimension, 
encoding sums of Feynman diagrams. This is made in Sec. |H In Sec. |5] we describe 
some tricks allowing to simplify the calculation of the Feynman diagrams. In the massless 
case these diagrams are calculated (exactly to all orders) in Sec. El In this section the 
sum of diagrams is reduced to integrals of total derivatives. From these integrals of total 
derivatives the exact NSVZ /3-function is derived exactly to all orders in Sec. [3 A 
similar investigation for the Pauli-Villars contributions is made in Sec. [HI The result is 
also given by integrals of total derivatives. The obtained expressions are verified by the 
explicit three-loop calculation in Sec. M The results are summarized in the Conclusion. 
Some technical details are presented in the Appendixes. 

2 = 1 supersymmetric electrodynamics and its reg- 
ularization by higher derivatives 



The massless = 1 supersymmetric electrodynamics is described by the following 
actionj^ 

S = ^Re j d^x d^e WaC'^Wb + ^Jd^x d^e (^0*6^^0 + 0*e-2^0) . (2) 

Here and are chiral matter superfields, and is a real scalar superfield, which contains 
the gauge field as a component. The superfield Wa is a supersymmetric analog of the 
gauge field strength. In the Abelian case it is defined by 

Wa = ^D^DaV. (3) 

[Da and Da are the right and left supersymmetric covariant derivatives, respectively.) 

In order to regularize the theory we add to the action a term with higher derivatives. 
Then the action can be written as 



Sreg = ^Re J d^x d^e WaC'^R + ^ J d^x d^e (^0*6^^0 + 0*e-2^0) , (4) 

where the function R satisfies the following conditions: 

R{0) = 1; R{oo) = oo. (5) 
For example, it is possible to choose 



^In our notation ri^i, = diag(l, —1, —1, —1); 9"^ = 0^0^°'; 6a and 6a denote the right and left components 
of 6, respectively. 
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In the Abelian case the superfield Wa is gauge invariant, so that action (jl]) is also gauge 
invariant. 

Quantization of model dl]) can be made by the standard way. For this purpose it is 
convenient to use the supergraph technique, described, for example, in textbooks [25| 126]. 
and to fix the gauge invariance by adding the following terms: 

-6i? / ivD'-D'R{^)v + VD'D^-R{^)v). (7) 

After adding such terms, a kinetic term for the superfield V will have the simplest form 

Sre, + S,i=^J d^x d^e Vd^R[^) V. (8) 

In the Abelian case, considered here, diagrams with ghost loops are absent. 

Adding the higher derivative term does not remove divergences in the one-loop dia- 
grams. In order to regularize them, it is necessary to insert Pauli-Villars determinants 
into the generating functional p51j. Therefore, the generating functional can be written 

as 

Z= I WZ^0D0 J](detPF(V^,M,))''exp(^(5,eg + 5gf + 5source)), (9) 

where 

^Sourse = j d'^x d^O J + { j d^x d^O {(pj + 4>~]) + h.c.) . (10) 

(It is necessary to substitute e in S'reg and S'gf by the bare coupling constant Cq.) The 
Pauli-Villars determinants are given by 



wher^ 



(^det PV(V,M)) ^ = J D<i> D<l> e'^P-^ , (11) 



Spy = ^J d^x d^e {^^*e'^^^ + <l>*e-2^<l>) + y d^x d^O M$$ + h.c.) . (12) 
The coefficients c/ satisfy conditions 

5^c, = l; J]c,M2 = 0. (13) 

Below we will assume, that M/ = a/A, where a/ are constants. Insertion of the Pauli- 
Villars determinants allows to cancel remaining divergences in all one-loop diagrams. 

The generating functional for connected Green functions and the effective action are 
defined by the standard way. 



^Note that the Pauh-Villars action differs from the one used in [21] because here the ratio of the 
coefficients in the kinetic term and in the mass term does not contain the factor Z . Using terminology 
of Ref. [32], one can say that here we will calculate the canonical coupling etc, while in Ref. [21] the 
holomorphic coupling ah have been calculated. Certainly, after the renormalization the effective action 
does not depend on the definitions. However, the definitions used here are much more convenient. 
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3 Generating functional 



In order to derive an exact /3-function we perform explicit summation of the corre- 
sponding Feynman diagrams. However, for the rigorous proof it is desirable to obtain 
some formal expressions encoding sums of these diagrams. For this purpose we first per- 
form the integration over the matter superfields in generating functional (|9]). This can be 
made, because the corresponding integral is Gaussian. It is easy to see that the result is 



Z = y^V^ JJ (^det PViy, Mi)y' det(*) det{*) 



X exp <; . / d'x ( ^,Vd'RidyA')V - j — * —f - J^*^f 

-0 



* 492^ 



(14) 



where 



^2V 



~2V 



(15) 



encode chains of propagators, connecting vertexes with the quantum gauge fieldH In order 
to obtain generating functional (JH]), we note that a solution of the motion equation for 
the chiral superfield 



8r 



(16) 



can be written as 



n2r)2 . _i r)2 



16^2 ^'J 2d^^ 16^2 * 2^2-^ ■ ^^'^^ 

Similar to derivation of Eq. f|T^ it is possible to perform integration over the Pauli- 
Villars fields. Introducing sources for the Pauli-Villars fields in expression flTTj) . we can 
write the Pauli-Villars determinants in a similar form: 



detPViV,M,j^^)-- 
where we use the following notation: 



1/2 

det(*) ) exp 



pv I 5 



A = 



-k is defined by 



4^2 



/ ^2 








\ 









^2 










f 








Z}2 





; 3 = 


J 


V 








D2 J 




\f / 



(19) 



I- I 







with 



'0 



VP, 



(20) 



■^This is a rigorous definition. In Ref. j29| * was defined only qualitatively using Feynman rules. That 
is why some equations below differ from the corresponding equations in Ref. |29] . 
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where 



V 



/ (e^^ - 1) 
fe^^ - 1) 



v 

























-2V 



-2V 



-1) 



-1) 



is a vertex contribution, and 
^ 



2 n2 



D'D 



MD' 



\ 



16(92 + M2) 

4(^2 + M2) 




16(92 + M2) 4(92 + M2) 








\ 





ML)2 





^2^2 



4(92 + M2) 

52^2 

16(92 + M2) 




4(92 + M2) 16(92 + M2) 

corresponds to propagators. 

The contribution of and can be also written in a form similar to flTH]) : 



Z= /" JJ (^detP\/(V,M/)y'(^det(*)y^^ 



X 



where M inside -k should be set to 0. 
It is also convenient to define 



(21) 



(22) 



(23) 



{h)a = [lo,0a]; (/l)a = [/o,^a]; (h) = {[lo, Oa], O'^h 

If, for examplfl Iq = {e^^ - l)D^D^/lQd'^, then 



D2n 

{Ii)a = {e -1)^^: 



ih)a = ie'^ - 1 



D D2 

892 



(24) 



(/2) = (e^^ - 1) 



492' 



{l2)ab = (e^^ - 1) 



492 ' 



(/a). = (e^^ - 1) 



292' 



(25) 



Let us also describe some properties of -k, which will be useful below: 
It is easy to see that 

-k iQ-k = —-k +-k^, 



(26) 



■*This expression can be used in the massless case, if only the field is present. In the general case Iq 
is defined by Eq. 
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and the following expansions take place: 

^=l + 5^(/or; ln* = 5Z-(/o)"; ^^ = 5^(n+l)(/o)^ etc. (27) 
Therefore, 



n 

n=l n=l n=0 



In^ „ = - )n = + 1) Wn.; )n = ^ 77 ^Wn; 



4^ _ (n + l)(n + 2)(n + 3) 



(^')n = ^ -^-^ (28) 

where the subscript n denotes the n-th term. 

4 Formal calculation of renormgroup function 



Expression ( fT4l) allows to obtain a simple formal expression for an anomalous dimen- 
sion. A Green function of the matter superfield is given by 



TTT^GidXy (29) 
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The corresponding inverse function, which by definition satisfies 

5^T Dl ( S'Y \ 1 



IS 



f-^G \d^)8l^ = -—— = (-^ * ^nsrO- (31) 



6(f)*J(l)J ' ' 53*J]^ \ 8^2 8^2 



ly 

(In order to derive the last equality, we have differentiated generating functional (fT4|) with 
respect to the sources.) The angular brackets denote the functional integration over the 
gauge superfield V . (The factors det(*) and det(*) should be certainly included.) From 
this equation we obtain 

DlDlG-\d')6%={*DlDl6i;). (32) 

This can be verified by applying to both sides the operator —DiDl/lQdy. We will use 
expression fl32|) later. 

It is also possible to construct an expression for a two-point Green function of the 
gauge superfield. For this purpose we use the equation 

r- ^ i,/ Va^W . i / ,«./,V.V.(.^|^| . (33) 
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which is derived in Appendix |X1 Here the symbol IPI means that in this expression it 
is necessary to keep only one-particle irreducible graphs, and Sj denotes the interaction. 
For convenience we denoted the argument of the effective action by the bold letter V. 
(This equation can be also easily obtained using the background field method, which is 
not used in this paper.) Here 



2V 



l)(/) + (/)*(e-2^-l)0j. (34) 
We substitute this expression into Eq. ( 133|) . taking into account the identity 

(^(^■^•» = I^(7^-7^)^- 

where the generating functional Z is given by Eq. ( IT^ . After simple calculations we 
obtaiijfl 



r(2) _ r,{2) 



where 



-5gf + (^-'- (Tr(Vg Jo*)) ' -2Tr(Vg Jo^Vg Jo*) -^Tr (V Vq*)) + (PV), (36) 



TrA^t.JS.A^„ (37) 

and tr denotes a usual matrix trace (if it is needed). (PV) denotes contributions of the 
Pauli-Villars fields, 



Q 



( 1 











\ 





1 
















-1 







^0 








-1 


/ 



Jn 
















\ 
































v 













/ 



P. 



{31 



Due to the supersymmetric Ward identity the two-point Green function of the gauge 
superfield can be presented in the following form: 



r 



(2) 



1 



16n J (27r)^ 



d^9Y{e, -p) ^^U^/2Y{9,p) d-\a, X,fi/p), 



(39) 



where a is a renormalized coupling constant, and 



S^Hi 



/2 



is a supersymmetric transversal projector. We will calculate the expression 



d 



c/lnA 



'\ao,A/p) - fto^) 



p=0 



daQ 
din A 



/3(ao) 



On 



(40) 



(41) 



(Here A and a are considered as independent variables.) From this equation it is evident 
that this expression is well defined. (Later we will demonstrate this in the three-loop 
approximation explicitly.) Note that here we implicitly use the higher derivative regu- 
larization, because it allows to perform differentiation with respect to In A and set the 
external momentum p to 0. 



^In order to find the contributions of (/)- and (?i-loops it is necessary to set M = 0. 
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5 Some useful tricks and summation of subdiagrams 



In order to calculate expression ( 1411) we consider 



d 



p(2) Q Q 



dlnAV^^ " "*^V p=o' ^^^^ 

Making calculations in the limit p — )■ 0, where p is the external momentum, is possible, 
because the corresponding integrals are well defined in this limit. The higher derivative 
regularization and the differentiation with respect to In A ensure that there are no IR 
divergences. This agrees with the results of Ref. [33] that the IR region does not affect 
to the /3-function. In order to obtain a transversal part of the two-loop Green function of 
the gauge superfield by the simplest way, we make the substitution 

V(x, 6) ^ ^"^a^^^6 = e\ (43) 

so that 

d^eY{x, e)d'^Ilii2^f{x, 9) -8. (44) 

(This is possible, because in the limit p — > the gauge superfield V does not depend on 
the coordinates .) In the momentum representation 

V(p,e) = j d^xV(x,^)e~'P""" ^ (27r)V(p)^^ (45) 
Thus, after substitution (H3|) we obtain 



d 



d 



rflnA 



-1 



p=0 



i27rY6\p) 



an 



p=0 



din A 



.(2) 
V 



s-s, 



gf 



p=o,v{x,e)=e'^ 



(46) 



We will try to reduce the sum of Feynman diagrams for the considered theory to 
integrals of total derivatives. In the coordinate representation such an integral can be 
written as 



Tr^[x'^, Something] j = 0. 



(47) 



In order to find a /3-function one should consider the massless theory. In the massless 
limit the fields and decouple. The Pauli-Villars contributions (for which this is not 
true) will be considered later. First, we will find a contribution of the field to the 
/3-function. The contribution of the field can be found similarly. We will take it into 
account in the end. 

In order to extract commutators (1471) . we consider diagrams containing a vertex to that 
only one external line (and no internal lines) is attached. We can add such a diagram to 
a diagram, in which the external line is shifted to the nearest vertex. Let us formulate 
this rigorously. In the massless case 
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,„^.« P^^; V^(e--l). (48) 

SO that Jo = Iq + P. As a consequence, 

*VJo = -^V(/o + P) = VP + -^f/oVP + V/o) = VP + *vfPVP + VP). (49) 

The expression V(PVP + VP) corresponds to a sum of subdiagrams presented below. 
Making a substitution V — t- 9°'9a0^9i, we obtain 



- (e -1)^-9 9J -^ + 9 0a^,+ie (7 \ e^^^ 19 (7 ). + 

(50) 

Only the first and the third terms give nontrivial contributions to the two-point func- 
tion of the gauge superfield, because they contain 9. Really, finally it is necessary to 
obtain 

d^9 9°'9a9^9b, 



while calculating a 0-part of a graph can not increase degrees of 9 or 9. Therefore, we 
should have 9°'9a from the beginning. 

6 Reducing the sum of diagrams to integrals of total 
derivatives 

6.1 One- loop approximation 

For the general renormalizable = 1 supersymmetric Yang-Mills theory, regularized 
by higher derivatives, the one-loop /3-function was calculated in |24]. The result is given 
by an integral of a total derivative and agrees with the exact NSVZ /3-function. Therefore, 
below we can make calculations starting from the two-loop approximation. 

6.2 External V-lines are attached to different loops of the mat- 
ter superfields 



Let us try to find a sum of Feynman diagrams exactly to all orders of the perturbation 
theory. We will start with diagrams in that the external lines are attached to different 
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loops of matter superfields. Let us consider a loop of matter superfields with n vertexes. 
This loop is proportional to 



Tr(#(7^)/^.(e^^ - 1)^^ * -9^9,9\e^^ - 1)^ * )^ = 

= '-Tr{^nY)c%{e^' - 1)^^ *^ -^^^.^-(e- - 1)^ *^ )^. (51) 

contains n — 1 vertexes, and one vertex corresponds to explicitly written (e^^ — 1).) 
After simple algebraic transformations this expression can be written as 

iTr( - e^ej" . (e- - 1)^ * -B-'e- . (e- - l)^ * (e- - * + 

* (e"' - 1)^^^^ * +e^e^9^e» terms) = 
= -Tr( - 29'9J''[9,, *] + tnr)c%[y;, *] + terms) , (52) 

n \ In 

where yl = — i9"'{piv)a9h is an antichiral coordinate, such that 

DaVl = 0- (53) 
Therefore, the considered loop with an external V-line can be written as 

1i(^-29'9,9''[9dM*]+i9\^^),''9d[y:M*] + 9\9\9\9^ terms). (54) 

Total derivatives in this expression give 0, because In* does not contain d^/d'^. Thus, the 
loop gives only 

9^,9\9\9^ terms. (55) 

The considered diagrams contain two loops of the matter superfields with an attached 
external V-line. Therefore, in order to calculate such diagrams it is necessary to multiply 
expressions ( l55l) . corresponding to each of these loops. As we explained above, in order 
to obtain a nontrivial result, it is necessary to have at least the second power of 9 and of 
9. Therefore, all these terms vanish after the multiplication and subsequent calculation 
of the diagram. Thus, the sum of all such diagrams is given by an integral of a total 
derivative and is equal to 0. 

6.3 External V-lines are attached to a single loop of the matter 
superfields 

Now let us consider a case in that both external y-lines are attached to a single loop 
of the matter superfields. For simplicity we will consider only a contribution of the field 
(f). (This means that the external lines are attached to the (p loop. Other loops can 
certainly contain 0-propagators.) The matter is that in the massless limit the fields 
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and (f) decouple. It is easy to see that a contribution of the field is exactly equal to the 
contribution of the field (f). 

In the massless case we should calculate the diagrams 



(56) 






where the vertexes are given by the corresponding terms in Eq. fl50|) . 

Let us shift ^-s to an arbitrary point of the loop, commuting them with matter prop- 
agators. This gives (the coefficients correspond to the expression fHB]) : contribution of 
is not taken into account r' 





2^2 



d 



TiieUe^^ -I 



64cilnA 



2(7n/^^Tr(ri(e 



* e 



2V 



+ (e2^ - 1 



c?lnA 



2V 



-2i- 



d 



-Tt(9 



2V 



c/lnA 

Jjdjj2 



Jjdjj2 
Jjdjj2 



' * (e^^ - 1 

2V .1) 



* e 



Q2 



* e 



2V 



^2V 



452 



892 452 



^2^2 ^ ^ 492 



^8^2 

jjd-Q2 

* 



* + e 



21/ 



1) 



492 



* (e' 



2V 



* e 



2y 



1) 



D2n 



* e 



2V 



1) 



jjd-[^2 



8^2 ' ' 8^2 ' ' 492 

We will start with the calculation of the following sum of diagrams: 



1 

+ 2 





(57) 



Using the identity 



5 



after simple algebraic transformations we obtain 



--^] = -2n'6\pE) = -2ti'i5\p) 



(5^ 



2i 



dlnA \ V 8 



(e2^ - 1)52d254(9„) + 



^In order to obtain the contribution oi a, (f) loop it is necessary to make a substitution e^^ 
and 



(59) 
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Terms proportional to the 5-function will be calculated in the next section. (So far we 
have not yet found all such terms.) We will see that they give a part of the /3- function 
proportional to the anomalous dimension. 
Now let us calculate the diagrams 




+ 



This sum can be written as 




(60) 



^Tr(^^((7'^)/[y;,(/0e*(/i)'*] +2(7^)/(/o)|* ((/i)c*(/i) 

+ {hy * (/i)e * ) - 2t(2{l2) * {hY * ih)d * +2(/2)/ * (/l)c * ihY * 



+4(/i)rf * {IiY * (/i)c * {hY * ))) + terms proportional to a ^-function. (61) 

(Terms proportional to a 5-function will be calculated later.) 

In order to present this expression as an integral of a total derivative we will use the 
identity 



Tr (^^ [{Yriy;, A]K B}[e,, c} + ii.n-^f' [Oa, B}[e,, C}[y;, A] 

-4z[r, [Oa, A}}[e\ BjlOb, C})) + cyclic perm, of A, B, C 



Tiie^^^Y' y*A%B}[ea,C} + {-lY-[9a,B}%,C}A 



+cyclic perm, of A, B, C, (62) 

which was proved in Ref. For the completeness we also present this proof in the 

Appendix [Bl Here A, B, and C are arbitrary differential operators, constructed from the 
supersymmetric covariant derivatives, which do not explicitly depend on 6, and Px is a 
Grassmannian parity of X. 

Qualitative arguments presented in Ref. [29| allow to suggest that expression (!6T|) for 
a diagram with n vertexes on the considered matter loop can be written in the following 
form: 



d 



dlnA \ 



Tr(e\YV\y:,{h)c*{hY 



6 



d 



n{n + l)(n + 2) d\nA 



Tile 



HYrK, *1 [Ob, *] [0a, lo] + {YYTa, *1 [0b, *] [y;, h] - ^0", [0a, A } [0\ *] [0b, h] 
HrY'[y;, *][0b, h][0a, A + W'[0a, *][0b, h][y;, *'] - 4z{r, [0^, *]}[0\ hWb, A 



+ {lT[yl, h][0b, A[0a, *] + Wya, h][0b, *] - 4z{r, [0a, hW, *'][0b, * 



a\ab 



31 r„ * 



36 ,31 



HrY'[y;, *][0b, *][0a, *] + {rY'[0a, *][0b, ^M, *] - 4z{r, [0a, *]}[0\ *m 



b, * 



(63) 
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In order to prove this, it is necessary to calculate all commutators and take into account 
Eq. (1251) . For example, if n = a + 6 + c + 3, then 



6 



A{*\B{*),C{*), + A{*),B{*%C{*) 



n{n + l){n + 2) 



(64) 



because 



^(a + 6 + c + 3)(a + 6 + c + 4)(a + 6 + c + 5) = ^(a + l)(a + 2)(a + 3) 

D O 

+ ^(6 + + 2)(6 + 3) + ^(c + l)(c + 2)(c + 3) + \{a + l)(a + 2)(6 + 1) 
o 2 

+ ^(6 + + 2)(a + 1) + i(a + l)(a + 2)(c + 1) + ^(c + l)(c + 2)(a + 1) 

+ ^(c + l)(c + 2)(6 + 1) + ^(6 +!)(& + 2)(c + 1) + (a + + l)(c + 1). (65) 

(A similar, but larger identity can be written for the term with four * in Eq. ( |6T]) .) 

Applying identity ( 162|) to (!63|) . we present expression ( 161]) as an integral of a total 
derivative: 



d 



dlnA 



Tr(^^(7^)/[y;,(/i).*(/i)'^*])^ 



n(n + 1)(?2 + 2) rflnA 



Tr(r(7 



X 



+terms, proportional to a (5-function. 
Calculating the commutators with 9 and ^ we obtain 



(66) 



d 



c/lnA 



2(7. 



Ii\ab 



Tr- 



n(n + l)(n + 2) (ilnA 



X{h)a + ( *^ (/l)a * + {h)a + * {h)a *^ ) * {h)b * h + *{h)b(^*'^ {h)a * + 

X(/l)a *^ + * {Il)a *^ ) + *{h)a * {Il)b *^ +/o ( *^ {h)b * + *^ {h)b *^ + * {h)b *^ ) 

X * (/l)a * +(/l)a( *^ {h)b * + *^ (/l)6 *^ + * {Il)b *^ ) * + *2 (/i)6 *2 (/i)a * + * {h)a 

X *^ (A)?) *^ ) + terms, proportional to a 5-function. (67) 

J / n 

Thus, the sum of all remaining diagrams is also given by an integral of a total deriva- 



tive. 
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In order to simplify the obtained expressions we derive an identity, which corresponds 
to shifting a loop momentum in an integral of a total derivative. For this purpose let us 
formally assume that y* and 6^ do not commute. Then due to the Jacobi identity 

[[d\y;],A] = [e\[y;,A]]-[y;,[e\A]]. (68) 

As earlier, we assume that A is a differential operator constructed from the supersym- 
metric covariant derivatives. As a consequence 

[y;,A] = -2^{J,r%[9,,A} + 0{9'). (69) 

Therefore, 



[[e\y;],A] = -2ti^,r[e\ea[e,,A}] + 2ti^^rX[0t, [0\a]} + 0(6') = o{e'). (ro) 

(Terms that do not contain 9'^ vanish after integration over d^O.) So, without using the 
relation [y*, 6^] = we obtained 

[[e\y;iA] = o{e% (71) 

Because the operation Tr includes the integration over dt^O, this means that it is possible 
to make cyclic permutations (Pa = Pb) 

Ti(e^[y;,AB]) = TT([9\y;]AB) = TT(A[9\y;]B) = (-1)^-Tr(e^[y;, (72) 



Actually this corresponds to shifts of the loop momentum in an integral of a total deriva- 
tive. Because the integrals are well defined, such shifts do not change the integral. 

Taking into account the possibility of making such cyclic permutations, one can sim- 
plify expression f l67|) . In Appendix [C] we prove that it can be written as 

a+b+2=n 

-|-terms proportional to a 5-function. (73) 
Collecting all results we obtain the following expression for a ^-contribution to f l46|) : 



■ ,D^D^d. ^ 2(6 + l)(7^)^'^,^ ■ 
?/;'2<e2'^-l)^^*- E ^ ^^^(/l)c(*)a(/l).(*). 

a+6+2=n 



dlnA 

-|-terms proportional to a 5-function. (74) 



The terms written explicitly are certainly equal to 0. The exact NSVZ /3-function is 
obtained from the terms proportional to a (5-function, which are calculated in the next 
section. 

Taking into account a possibility of making cyclic permutations in the expression which 
is commuted with |/*, ( !7i|l can be rewritten as 
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d 1 
din An 



Tt(9 



2V 



D'D'd, 
' 16^4 



*-2(7'^)^^*(/0e*(/l)d* 



+terms proportional to a 5-function. 



It is easy to see that this expression can be presented in the form 
d i 



-Trie 



{y^y, [yl, *] ) + terms proportional to a ^-function. 



din An 

Thus, (taking into account contribution of 0-loops) finally we obtain 
. d 



dlnA 



-Trie 



(yr^KMi**)] 



+ terms proportional to a (5-function. 



(75) 



(76) 



(77) 



7 Derivation of the NSVZ ^-function 



In the previous section it was found that all integrals giving the /3-function are integrals 
of total derivatives. However, they are not equal to 0, because 



{f{oo) = due to the higher derivative regularization.) This is equivalent to taking into 
account terms with a 5-function. Really, let us rewrite this equality as follows: 



oo 


= -J d'qS\q)f = -fiO). (79) 

Thus, we see that total derivatives with respect to q^ are equivalent to total derivatives 
with respect to q^ plus terms proportional to 5^(g). In the approach described in the 
previous section 5- functions appear, if is commuted with 9^/(9^ as in Eq. ( 158|) . In this 
section we calculate all such terms. 

Qualitatively, the (5-function allows to perform an integration over a momentum of the 
considered matter loop. This corresponds to cutting the loop, which gives diagrams for 
the two-point Green function of the matter superfield [22]. For example. 




Let us derive this by a rigorous method. First we consider expression (1591) . Omitting 
total derivatives (corresponding to d/dq^) we obtain 
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d 



j d\d^5l^{'^e^ * (e- - 1)( - ^)D^DH\d)5l^ 



d 



cilnA 



4 n2 n2x4/QN x8 



a: J/ 



(80) 



Taking into account that 



dS 



s\d)s\x -y)= -^5\q)e-"i-^^''-y''^ 



(2vr) 



and calculating 6'-integrals, we obtain 



.^TJ!lleMe'''-l)D'D'S\d)) = f,,J, 
rflnA \ 4 ^ ' ^ V (27r)4dlnA 



d'xd'y(*DlDlS% 



12) 



Using Eq. (132!) for the function G ^ this can be presented in the form 



TT 



d 



(27r)4rflnA 



47r2 

~ (27r)4rflnA 



d^xd''yG-'6\x-y) 



An' 



d 



dlnA 



p=0 



^3) 



This expression is not well defined. Thus, it is written formally. However, we will see that 
after adding the other contributions a well defined result is obtained. 

5-functions are also present in expression (1731) (or fl6ip ). if the matter loop contains 
coinciding momentums. Really, taking into account that 



(e^^ - 1 



D^Dr 



„2V 



1) 



4) 



d^/d^ appears due to the following identities (It is assumed that momentums in Ji and 
7i coincide): 



(/i),-(7i),^ 



852 8^2 



[h)d ■ [h)c 



^((l + 75)7'^)c.^/^^/^^; 



^5) 
^6) 



Two (or more) momentums coincide, if two cuts of the matter loop make a diagram 
disconnected. An example of such a diagram is 
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For analyzing such diagrams we will use the identities 



Di^D"^ D'^D'^ DhD^ , D^D^ D'^D„ D^D„ 



(T'\ T ^ U^_U_ _ _U^U^ , , U U U Ug 

UJfe • ^ ■ 16^2 - - 8^2 ; ■ UJa ^ ^^g^2 ' 8^2 " " 8^2 



^o-^^3^-^^--^^- (87) 

Let us assume that there are p coinciding momentums q in the considered matter loop 
(to that the external lines are attached). Then the corresponding diagram contributing to 
the (connected) two-point Green function of the matter superfield is not IPI and consists 
of p parts, connected by a single line of the matter superfield. (If there are several groups 
of coinciding momentums, each group should be considered separately.) 

Let us assume that the parts of such a diagram contains q (z = 1, . . . ,p) vertexes on 
the matter line (to that the external lines are attached) III We will denote expressions for 
these parts by Gi, G2, • • •, Gp- Due to identities flHTj) the following variants are possible: 



a + 1 = Ci; 6 + 1 = C2 + C3 + . . . + Cp; 
a + 1 = C2; 6 + 1 = ci + C3 + . . . + Cp; 

a + l = Cp] 6 + 1 = ci + C2 + . . . + Cp_i, 
where = ci + C2 + . . . + Cp, (88) 

because terms with (/i)-s give a nontrivial result only if there are no lines with the 
momentum q between (/i)a and {Ii)b- (We assume that momentums in (Ji) and (Ji) are 
equal to q.) However, any number of such hues can be between {Ii)h and (/i)a- According 
to the results of the previous section, it is necessary to calculate (and subtract) a singular 
part of the expression 



a+b+2=n 

A singular part of the first term has been already found and is given by (183|) (with the 
opposite sign). For a diagram that contains a sequence of subdiagrams Gi, G2, • • it 
can be written as 

^For the diagram, presented above, p = 2, ci = 2, and C2 = 4. 
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d 



dlnA 



Tt 9 



16(9^ 



Singular part 



G1G2 ■ ■ ■ Gp^ 



\p-i 



P 



71^ d 

4 cilnA 



i:i[e^GiG2...GpD'D'8\q)]. (90) 



The factor p is present, because there are p variants by which D D d^/lQd can be placed 
between Gj. 

In order to calculate a singular part of the second term, we consider 



2(7 



li\cd 



d 



dlnA ^ — ' n 

a+b+2=n 



- Singular part 



(91) 



Using a possibility of making cyclic permutations inside the commutator and identity 
( l86l) . this expression can be written as 



-1 



\P-2 



d 



-Trie 



y;, ^,G,G2 . . . G,DW'D, 



X 



rflnA 

Ci + . . . + Cp-i Ci + . . . + Cp-2 + Cp 



Ci + C2 + . . . 



+ 



Ci + C2 + . . . c„ 

(-ir(p-i) 



+ . . . + 



Singular part/ 

C2 + ■ ■ ■ + Cp 
C1 + C2 + ...Cr, 



TT^ d 

TrfhTA 



TT[e^GiG2...GpD'D'6\q)), (92) 



because 



D'^D^DaS^d) = D^DH\d). (93) 
Therefore, a ratio of the coefficients in the first and the second terms of Eq. (l89il is 

p — 1 



p 

Taking a sum of both contributions we see that a coefficient is proportional to 

, 1 1 



(94) 



(95) 



p p 

The coefficient 1 corresponds to the expansion of G~^ (see Eq. ( |83l) ). Therefore, taking 
into account the expansions 



ln(l-a:) = -V-; = E 



(96) 



p=l p=0 

we see that singular parts of the commutators give the contribution to Eq. (H6 



^T^'b\p) 



d\nG 
=0 c? In A 



-47r'r(p) 



(97) 
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Here the expression 7(00) is well defined, unlike the corresponding expression in Eq. fl5^ . 
Thus, after taking into account all contributions, the well defined result is obtained. 

Diagrams with a loop of 0-fields (to that external lines are attached) give exactly the 
same result. Therefore, due to Eq. (146!) a /3-f unction is given by the sum of the one-loop 
contribution a^/vr and 

A/3 = -— 7(«). (98) 

TT 

Thus, we obtain the exact NSVZ /3-function 

/3(a) = -(l-7(«)). (99) 



8 Pauli— Villars contributions 



Previous calculation was formal, because so far we did not take into account contri- 
butions of the Pauli- Villars fields. However, these contributions can be considered in a 
similar way. 

8.1 Summation of subdiagrams 

Let us start with the summation of subdiagrams. For the Pauli- Villars fields there 
are four different types of subdiagrams. Below we will calculate elements 22, 23, 14 
and 11 of the corresponding matrix. After the substitution V ^ 6"^ and some algebraic 
transformations (omitting for simplicity expressions for the left vertexes) they can be 
presented in the following form: 



+ 



+ 



(92 + M2)2 " 4(92 + M2) 



16(^2 + M2) J 



terms without Q (100) 
+ terms without Q. 



' 16(^2 + M2) 



+ 



+ 



+ 



MDh 



2(92 + M2)2 92 + M2 

ML>2 . 

y;. ...o . .... -2W 



4(92 + M2) J 



terms without 9 



(101) 



b, 



4(92 + M2) J 



terms without 9. 
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e 



-"2(92 + 1/2)2 " 4(92 + M2)2 

ML)2 



^l■> 



4(^2 + M2) J 



9^ 



+ terms without 6 



(102) 



4(^2 + M2) J 
-^''-TTT^ + teims without 6. 

4(^2 + M2)2 



+ 



+ 



8(^2 + M2)2 



anb DaDb 



,an nb 



92 + M2 



+ 9^9 J 



4(^2 + M2) 



(103) 



,5 D,DW^ 



9" 



-l9\l^,)a% 



16(92 + M2)2 92 + M2 

^^^^ 1 ^^a...r. D-D- 



+ terms without 6' 



16((92 + M2) 



^6, 



-9' 



DkD'D 



2 n2 



16(92 + M2) 

.-6 ^6 



16(92 + M2)2 92 + M2 



terms without 9. 



The other matrix elements are calculated similarly. The whole matrix corresponding to 
sums of subdiagrams fllOOj) — fll03p is written as 



Q (^"V X 
/ 



+ 














MDaD 



\ 

2 



4(92 + M2)2 92 + M2 



MDaD- i{r)a''D,D-d^ 



4(92 + M2)2 









D. 



4(92 + M2)2 4(92 + M2)2 92 + M2 







+t9''{Y)a%[y* Qh] - 29''9j%, QIo] + terms without 9 



where 



(104) 



/ -1 \ 

10 

10 

\ 0-1/ 



(105) 
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satisfies the identities 



[Q,Io]=0; [Q,*] = 0; = I. (106) 

Jo is defined by Eq. fl20l) . and V (corresponding to the left vertexes) is given by fl2T]) . 
Terms without 6 do not contribute to diagrams for the /3-function. (The integral over d'^O 
is nontrivial only if a diagram contains 6^ and, in particular, ^^.) 

8.2 External lines are attached to different matter loops 

This case is very similar to the massless one. A loop of matter superfields to that an 
external line is attached is now proportional to (for simplicity we omit — ^ c/) 



TT(^ie%Y)c%[y;, Qio] ^ -2e'e,e^[0d, qio] ^ terms)^ = 

= -Tih9%Y)/9,[y;, Qio] -29'9,9''[9,, Q/q] +9^ terms) . (107) 

n \ ^ In 

After some simple (but nontrivial) algebraic transformations this expression can be rewrit- 
ten as 

^Tr(Q ( - 29'^9,9'' ^ \9^, Jq] * \%Hl")cQ<i * b;, ^o] ^ ) + 9\9\9\9° terms)^ = 
= ^Tt(^q(^ - 29'9j''[9d,*] + ^9%J'')/9,[y;,*]) + 9\9\9\9<' terms)^. (108) 
Therefore, as earlier, a matter loop is given by 

TT(^Q(^-29'9j''[9dMi.]+t9%Y)c%[y:Mi<]) +9\9\9\9' terms). (109) 

As in the massless case, multiplying expressions for two such loops we obtain that all 
diagrams in that external lines are attached to different matter loops are given by integrals 
of total derivatives. All these integrals are evidently equal to 0. 

8.3 External lines are attached to a single matter loop 

Calculation of such diagrams in the massive case has some differences from the massless 
case. We construct Feynman rules in the massive case using Eq. (I104p . Note that they 
are different from the corresponding rules in the massless case, because the expression 

^9^ma%[y;Jo] (110) 

in the massless case gives 

^ni')aX [y;. ^] = ^ni^aX^^^ - e^o\^ (m) 
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and contains terms, proportional to O'^d^Oi,. 

Also in the massive case it is necessary to take into account the effective diagrams 




(112) 



The first diagram corresponds to terms, proportional to 9°" in Eq. (11041) . The second one 
contains a sum of subdiagrams with two adjacent external lines. These subdiagrams are 
presented in Appendix 

Let us now write down the results for all diagrams (again, omitting — ^ c/ for sim- 

phcity): 

The diagrams contributing in the massless case are calculated similarly. Taking into 
account identities (I106p . the result can be written as 




where 

^^^^^^rfi^^K^'t^'^^i^t^'^^]^)' 

Ex, = 2(7'^)/^Tr(^^ ([|/;, (Ji) J * (hY * h] ^ {h)c ^ {hY * ) ); 

= -2^^i^Tr(^4 ( - (Ji), ^ {h) ^ ChY * * {hY * * ih)c * {hY * 

+2(/i)d ^ {hY * (/i)c ^ {IiY ^ ) ) ■ (113) 
Taking into account that 

* [OaJo]* = [0",*]; Tr(Ar,5]) = Tr({A^n5) (114) 
and using Eq. (I106p . we also obtain 




(115) 
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4 ij-ff^Y'D^Dbd^ 
92 + M2 (92 + M2)2 




V 












ML)2 









\ 

ML)2 



(92 + M2)2 





4 ti-ff^y'DgD.d^ 

(^2 + M2)2 92 + M2 (92 + M2)2 







The last diagram is calculated using expressions for the subdiagrams with two adjacent 
external lines, presented in Appendix The result is 




d ^ , 4M2 ^ 



din A 



+i 



d 



dlnA 



(a2 + M2)2- 

/ 




Tr(r *V 





^ MD2 



\ 

^ MD2 



(92 + M2)2 












V (92 + M2)2 

Now let us find a sum of these diagrams. Similar to the massless case 



(116) 



Exi 



i d 



2dlnA 



din A 



/ 



d 



—l- 



dlnA 



Tr( e^V 



where we take into account that 



(92 + M2)2 92 + M2 
















(92 + M2)2 ° 




4 

(^2 + M2)2 ~ ^2 +M2 / 







(117) 



16(^2 + M2) 



2(^2 + M2)3 + (92 + M2)2 ^ 92 + M2 



9U - 



+ 



2(92 + M2)3 (92 + M2)^ 



92 + M2 



:ii8) 



Singularities, giving (5-functions, are certainly absent in the massive case. 
Summing ( IllSp . (]116p . and ( 1117^ we obtain 



^^Tr(^^^[y;,b;,/oM) = 0. (119) 
Moreover, exactly as in the massless case for diagrams with n vertexes on the matter 



loop 
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d 



n + + 2) c/lnA 



0. 



(120) 



(It is evident that there are no terms proportional to a ^-function in this case.) Taking into 
account a possibihty of making cychc permutations (172|) . we obtain that the contribution 
of the Pauh-Villars fields is given by the following integral of a total derivative (for 
diagrams with n vertexes on the matter loop to that external lines are attached): 



d 



dlnA \ 



a+b+2=n 



0. 



1211 



As earlier, the left hand side can be rewritten as 

i d 



2d\nA 



-Tt(9 



0. 



(122) 



The final result for the sum of diagrams in that the external V-lines are attached to 
a single loop of matter superfields is 



d 



din A 



-Tr( e' 



(y'^)*,ln(*) + ln(T)-l5^c,ln(^,) 



— terms with a (5-function, 



(123) 

where -kj means that it is necessary to use the mass Mj in the definition of -k. The 
factorization of integrands into double total derivatives, which follows from this equation, 
agrees with the arguments presented in Ref. [22] . 

From Eq. f ll23p it is possible to explicitly construct an integral of a total derivative, 
which is equal to the sum of Feynman diagrams. In the three-loop approximation this is 
made in the following section. 



9 Three-loop verification 



In order to verify the expressions, obtained in the previous sections, we compare them 
with the explicit three-loop calculation, made in Ref. [21] by a different method. The 
result can be written as [29] 



d 



din a 



d ^(ao, A/p) - ftg 



-1 



p=0 



16n{Ai + A2 + A3 



(124) 



^This result was also presented in Ref. |34j . but some sign in A3 were written incorrectly. 
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A 



ci 



d'^q Id d 



(27r)4g2 dq^dlnA 
d^q I d d r d^k 1 



(^In(g2 + M2) + 



Ml 



g2 + Mf 



(125) 



{2TiYq^dq^d\nK J {2ttY Rl\{k + q) 



i y 



(g2 + M2)2 



1 



X 



{{k + g)2 + Mj 
I 



dH 



2e2 



(27r)4 t\k + tf 



2e2 



+ 



(27r)4 {t^ + M^){{k + tf + Mj) 
d^q I d d f d^k dH 4e^ 



(126) 



(27r)4 g2 dq^ d\nk J (27r)4 (27r)4 Pi?^ PRi \ {q + k) 



2{q + iy 



2P 



(g + A; + /)2(g + /)2 (g + A; + /)2 
1 

+ 



^ ""'Xq^ + M2)2 ((g + kY + Ml) 
2A;2 



2((g + /)2 + M2) ((g + A; + /)2 + M2)((g + /)2 + M2) 
1 



+ 



2M| 



((g + A; + /)2 + M|) ((g + A;)2 + M2)((g + A; + /)2 + M|; 



+ 



2M2 



+ 



2M2 



(g2 + M2)((g + /)2 + M2) ((g + /)2 + M2)((g + A; + /)2 + M|) 



(127) 



where Rk = R{k'^/A'^). Here Ai is a one-loop result. A2 is a sum of two-loop diagrams, 
three-loop diagrams with two loops of the matter superfields, and diagrams with insertions 
of counterterms arising from renormalization of the coupling constant. A3 is a sum of 
three-loop diagrams with a single loop of the matter superfields. 
The anomalous dimension can be written as 



7(ao) 



dlnZ 



d 



dlnA dlnA 



(^ln(ZG) -InZ) 



dlnG 



q=0 



dlnA 



(128) 



5=0 



where A and the renormalized coupling constant a are considered as independent variables. 
The two-point Green function of the matter superfield in the two-loop approximation is 
given by the following integrals [211 



InG 



d'^k 



2el 



X 



{2TiYk^Rk{k + qY 



2el 



■ 1 r dH 2el 

, 'TJ ^^¥{k + tY 



dH 
(2^ 



(t2 + M2)((A; + t)2 + M2) 
4 



d^k dH 



'0 



+ 



{27tY {27tY k^kPRi 
8A;2 - 4g2 



(g + fc)2(g + /)^ 



(q + kYiq + k + lY (q + kYiq + k + lYiq + lY 



(129) 



Therefore, taking into account one-loop renormalization of the coupling constant, we 
obtain 
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d^k d 



{27iydlnAk^Rl 
dH 2e2 



A 



dH 



2^ 



(27r)4 (i2 + M|)((A; + t)2 + M2; 



(27r)4 t2(A; + t)2 



+ 



(27r)4(27r)4 d\nAk^Rkl^Ri{k + l)'^' 

(130) 



This expression is finite both in the UV and IR regions. The UV finiteness is ensured by 
the regularization. The integral is IR finite due to the differentiation with respect to In A, 
which should be performed before the integration. 

Taking the integrals of total derivatives in Eqs. fll25p — f ll27p using the identity 



dS 1 d 



■(/(g2 = oo)-/(g2 = 0) 



we obtain 



As a consequence, 



/3(ao) 



1 / _ d\nG 



an 



n 



dlnA 



q=0 



+ 0{al) 



/3(a) 



a 



1 -7(a)) +0{a' 



(131) 



(132) 



(133) 



In order to verify Eq. f ll30p . it was compared with the result of the calculation made 
with the dimensional reduction. (Such calculation can be made using Eq. (11290 .) Using 
the standard technique of the dimensional reduction one obtain 



a a" ^. n 



TT 



(134) 



7DREDiaj 

This result up to notation^ agrees with the calculation made in [5]. 

From the other side, it is possible to calculate the two-point Green function of the 
gauge superfield using Eqs. (l89ll and fll2ip . Taking into account that a contribution of 
diagrams with a 0-loop is equal to a contribution of diagrams with a (j)-\oop, we obtain 



d 



d\nA \ 



16^4 



= -327r 
-2e 

+647r 



d'^q 1 d f d'^k d 



(27r)4g2rfg2 J (271^ d\nA k^Rl{k + q)^ 
dH 1 /■ dH 

I 

d^q 1 d f d'^k dH d 



+ terms with a 5-function 



A 

Rk[l + —\n 



(27r)4 t\k + ty 



(27r)4 {t^ + M^){{k + tY + Mf) 



(27r)4g2t/g2 J {2'nY{2TiY d\YiAk'^RkPRi\{q + kf{q + k + lf 

{2q + k + lf 



(135) 



{q + kf{q + iy{q + k + iy 



^In order to obtain the results of Ref. f5| it is necessary to set a = g^/47r, 7(a) = 27(g), /3(a) 
<7^(ff)/27r. 
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The last term in Eq. fl5^ is 



a+b+2=n 



_l_^(,.)-Tv(.'[,;.(A).(«).(A).(.) 



+terms with a 5-function 



= 327r 
+ 167r 



rf^g 1 /• d^k dH d 



{27t)^ q^dq^ J (27r)4(27r)4 rflnAPi^fcF^, (g + A;)2(g + /)2 



dS d^k dH d 



d 



{2q + k + iy 



(27r)4 (27r)4 (27r)4 din A k'^R^PRi ^q^' q^{q + kf^q + iy{q + k + 1} 
Contributions of diagrams with a Pauli-Villars loop can be calculated similarly: 



r.(136) 



d 



2^ ^rflnA 
J 

= 327r^cj j 



Tt(9 



2 A 



e , 
+ — ^In 



47r2 /i 



rf-^g 1 


d r 


(27r)4g2 


dq^J 


/• dH 


2e2 



r 



(27r)4g2t/g2 y (27r)4rflnAA;2/22 (g2 + JVf2)2((A; + g)2 + M2) 

, X . I 2e2 

(27r)4 t2(A; + t): 



1 



Cl 



(27r)4 (t2 + M2)((fc + t)2 + M2) 



.647ry c, / 



d^q 1 d f d^k dH d 



X 



X 



4 g2 dq^ J (27r)4 (27r)4 dlnAk'^RkPRi (g2 + m2)2 
1 (2g + A; + If + 2Mf 



i{q + k + ly + M2)((g + A:)2 + Mf) ((g + A;)2 + M2)((g + If + M^) 
1 



+ 



2M2 



((g + A; + 0' + M2) (g2 + ^^((g + /c)2 + M2)((g + /)2 + M^) 
2M| 



((g + fc)2 + M2)2((g + k + lf + Mj) 

d'^q 1 d f d'^k dH d 



-327rJ] 



Cl 



(27r)4 g2 rfg2 j (27r)4 (27r)4 dhiKk'^RkPRi (g2 + M2)2((g + A;)2 + M^) 



d^g (i^/c (i"^/ (i 



d 



(27r)4 (27r)4 (27r)4 dlnAk^kPRi dqf" (g2 + Mf) 



(2g + A; + If 



((g + A;)2 + M2)((g + If + M2)((g + k + If + M]) ' 



(137) 



E .,^/t^ 7T^Tr(g^ Z/:, (/l)e(^)a(/l).(^) 



c/lnA 



0. 



/ a+b+2=n 

Summing all these contributions with the one- loop result, we obtain 



(138) 



^%^ = 167r(Ai + A2 + A3). 



On 



(139) 
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Thus, the general results, presented above, agree with the explicit three-loop calculations, 
made by a different method. 

Also it is possible to verify expression fll23p . After a calculation of Feynman diagrams 
we obtained 



tr{| ln(.^^ 



-2e' 



dH 



(27r)4 (27r)4 fc2i?2(g2 + M2)((g + kf + M^) 
1 ^ f dH 1 



.2 



^In 



(27r)4 f^i^k + ty 



E 



cj 



d^q d^k dH 



(27r)4 {t^ + M^j){{k + tY + M}) 
q^-M^ 



(27r)4 (27r)4 (27r)4 k^R^PRi \ {q^ + M^iiq + kf + M^){{q + If + M2) 
(g + kf - M2 4 



+ 



(g2 + M2)((g + A;)2 + M2)2((g + / + A;)2 + M2) 3(g2 + M2)((g + A;)2 + M2) 
1 8 
^ ((g + /)2 + M2) ^ 3(g2 + M2)((g + A;)2 + M^){{q + k + lf + M^) 

(2g + A; + /)2 + 2M2 
~(g2 + M2)((g + fc)2 + M2)((g + /)2 + M^){{q + k + If + M^) 



(140) 



where tr is a usual matrix trace, which (unlike Tr) does not contain J d^x. As a conse- 
quence, after some simple algebra we obtain 



/^M ^ >^ f dS d d In(g2 + M2) ^ 
27r^^ — r/ cj I — — —Tr~~ n \- Att 



On 



d d 



rflnA 
1 



dq^ dqf^ \q'^{k + qf 



(27r)4 dq^ dq^ q^ 
1 



d'^q d'^k 



ci 



-2e' 

-f47r 
4 



(27r)4 t2(fc + t)s 



(g2 + M2)((fc + g)2 + M2) 
rf^t 1 



rflnA J (27r)4 (27r)4 Fi?; 
2 A 



Cj 



d\ d^k dH 



(27r)4(t2 + M2)((A; + t)2 + M2) 

4 9 a r/ 2k^ 



dlnAj {2^lY{2^lY{2^iYk^RkPRl^q^'^q^\\ q^{q + kf{q + iy{q + k + lf 
2 \ ^ f 2{e + Mf) 



q\q + ky{q + iy 
1 



+ 



(g2 + M2)((g + ky + M^){{q + If + M^) 
2 1 



((g + A; + /)2 + M2) [q^ + Mj){{q + kf + Mj){{q + If + M^) {q' + Mjf 



X 



((g + A;)2 + M2)((g + /)2 + M2) 
Then, using the equation 

d (q^ 



integrals of 5- singularities. 



(141) 



{2t:Y \dq^\q 



f{q) - 2n'5\q)f{q) 



dS 1 df 
(27r)Vd^' 



(142) 
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we again obtain 



1671(^1 + ^2 + ^3). 



(143) 



10 Conclusion 

In this paper we have proved that the /9-function in iV = 1 supcrsymmetric electrody- 
namics, regularized by higher derivatives, is given by integrals of total derivatives in each 
order of the perturbation theory. In particular, explicit expressions for all total derivatives 
are obtained. Having compared them with the explicit three-loop calculation (made by a 
different method) we obtained the complete agrement. 

Factorization of integrands into total derivatives is the origin of the exact NSVZ /3- 
function, because one of the loop integrals can be taken. After this the /3-function in n-th 
loop is related with the anomalous dimension in n — 1-th loop. This was also proved in 
this paper by exphcit summation of Feynman diagrams. 
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A Derivation of expression for the two-point function 

Splitting the classical action S into a quadratic part and an interaction Si, it is possible 
to present the generating functional Z in the form 



where the Pauli-Villars determinants are omitted for simplicity. Differentiating this gen- 
erating functional with respect to the source J we obtain 



Z — exp 



■'^'Zo, (144) 



5Z_ _ 

S^Z 



( lie 1 1 \ 



exp(i5'/)(-2ieo);^Jx^o = -2iel——J^Z - 2iel 





dlR5vS~i^'~i'^3 




(145) 



5Jx5Jx 



As a consequence 




1 5Z _ 
Z 5Jx 
1 6^Z 




1 6S1 1 6S1 
dlR5VxdlRWy 




(146) 
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(In the last equation we set J = 0.) Therefore, if we denote an argument of the effective 
action by V, the two-point Green function of the gauge superfield will satisfy 



5^T .SHnZ 



i- 



5V.x5V„ / 6JxSJy 6JxSJy 



- r4.eV^^^^\ + 4e^^^/^\^ (147) 



The inverse matrix is evidently given by 



'^^r 1 , /_.SSi6Si , 6^ Si 



^ ^d'RSly+ U—^—^ + ——^) , (148) 
6VJ\y 2el y \ 5Vx5Vy 5VJVy/im 

where the symbol IPI means that it is necessary to keep only one-particle irreducible 
graphs in this expression. Thus, a part of the effective action quadratic in the gauge 
superfield can be written as 

= A /" d^xYd'RY + \[ d'xdSYxVyU^4^^4^ + -^^) . (149) 
^ 4e2y 2 J " y\ 5Vx6Vy 6VJVy/im ^ ^ 

The Pauli-Villars determinants can be considered similarly. In this case Sj contains the 
Pauli-Villars fields, and ^ cj should be also included. 



B Prove of identity (1621) 



Let us consider 



X ^ iT[e'{{^^r'[y;,A][e,,B}%,c} + {^n'^\-iY-[ea,B}[e,,c}[y;,A] 

-4i[r, [Oa, A}}\e\ C}) } + cyclic perm, of A, B, C, (150) 

where A, B, and C are differential operators, containing supersymmetric covariant deriva- 
tives. It is important that they do not explicitly depend on 6. Certainly, we assume that 

PA + PB + Pc = 0{mod2). (151) 

Using the identities 

Tr([y;,A]i?) =-Tr(A[y;,i?]); 

{-lf-[[A,B},C} + {-inC,A},B} + {-lfc[[B,C},A} = 0, (152) 

we obtain 
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X = {Yr'Tr[e'' ( [y;, A[e,, c} + B}[e,, c}a\ - A[e,, [y;, b]} 

X [9a, C} - A[9,, B}[da, [y;, C]} - {-if^K [y;, B]}[e,, C}A - B} 
[9,, [y;, C]}A^ } - 4zTr A}}[e\ B}[e,, C}) + cyclic perm, of A, B, C. 

(153) 



The similar operation is repeated for ^-s in double commutators: 



X = {Yr'T:t[e'[[y;,A[e,,B}%,c} + i-if-[da,B}[e,,c}A\ + i-if-[9,,A} 

X [y;, B][9a, C} + {-lf-A[y;, B][9,, [9^, C}} - {-lf^[9a, A}[9,, B}[y;, C] 
-(-lf-A[9a, [9„ B}}[y;, C] - [y;, B][9„ C}[9a, A} + {-lf^[y;, B][9a, [9,, C}}A 
-{-lf''[9,, [9a, B}}[y;, C]A + {-lf-[9a, B}[y;, C][9„ A})) 

-4iTr(^^^[^", [9a, A}}[9\ B}[9b, C}} + cyclic perm, of A, B, C. (154) 
In addition to 9^, the commutators with give one more degree of 9: 

[y;, A] = -2z{Yr%[9t, A} + 0{9'). (155) 

That is why it is necessary to be careful commuting 9'^ with A, B and C. For example, 
taking into account that all expressions containing 9 in less than fourth power vanish, we 
obtain 



[-lY-[^^Y^T,9'[9,,A}[y;,B]\9a,C} = {-lY-^\rT'T^^^^^^ 
-29%9''[9d, [9a, C}}[9,, A}{-2i){^^y^9e[9f, B}) 

= l^[{-lY-'^^9\Yr\0a,C}[9,,A}[yl,B]+Ai9^[9\ [9a,C}}[9\A}[9,,B}). (156) 
Similarly one can derive the following identities: 

= Tr ( - 9\^rr\y;, C][9,, A}[9a, B} + m'[9\ [9a, B}}[9\ C}[9,, A}) ; 
{-l)^B^\YrTr(9'A[9a, [h.B}}[yl,C\) (158) 
= 'Ii[{-iy'-^'e\rr\Oa, [9,, B}}[y;, C]A - 4z(-l)^-^^[r, A}[9a, [9\ B}}[9,, C}) ; 
{-lf-{YrT^^{e'A[y;,B][9,, [9a,C}}) (159) 
= {-lY^l^(9\^^Y'[yl,B\[9,, [9a,C}}A + m'[9^A}[9,,B}[9\ [9a,C}}). 
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Using these equations, after some algebraic transformations X can be rewritten as 



X = TT{e\rr''(^\y;,A[e,,B}[e,,c} + {-if^[ea,B}[e,,c}A\ - 2[y;, A][e,, b} 

+cyclic perm, of A, B, C. (160) 
Comparing this expression with the definition of X, we obtain 



x = -2x+ (TT{e\^,r' [y;, A[e,, c} + B}[e,, c}a] ) 

+cychc perm, oi A, B , C \ . (161) 



Therefore, 



X = ^Tr (^^(7^)^'' [y;, A[9,, B}[9,, C} 
+cychc perm, of A, B, C. 

This completes the proof. 



-if^[e,,B}[ef„c}A 



(162) 



C Simplification of expression (l67l) 



Let us calculate 



-Tr( e 



, ^o^-\^ k,*'(/l).*(/l)c+(*'(/l)c* + *'(/l)c*' + *(/l)c*') 

X (/i)d * Jo + *(/l)d( *^ (/l)c * + *^ {h)c *^ + * {h)c *^ ) + *{h)c * {Ii)d *^ 
+/o( *^ ih)d * + {h)d + * ih)d *^ ) * {h)c * +(/i)c( *^ ih)d * + {h)d 



+ * (/i)d *M * + (/l)d (/l)c * + * (/i)c (/i)d 



(163) 



for a diagram containing n vertexes on the matter loop to that external lines are attached. 
Due to a possibility of making cyclic permutations ( I72p 



iYr^r{[y;,{hU*)aih)d{*h 

Therefore, using the identity 



-(7^)^'Tr( y;,{{lM*WM*)a 



(164) 



(165) 



and Eq. fl25]) . the considered expression can be written as 
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J2 Cahr'Tr{\y;,{lM*)a{IlU*)b 



(166) 



a+6+2=n 



(Two vertexes correspond to (/i)c and {Ii)d-) In order to find the coefficients c^, it is 
necessary to calculate all (*^)a and (*^)b, using Eq. fl28|) . The result is proportional to 



-ha + l)(a + 2)(a + 3) - ^(6 + l){b + 2){a + 1) - hb + l)(a + l)(a + 2) - ^(a + 1) 
6 2 z 6 

x(a + 2)(a + 3) + ^(6 + 1)(6 + 2){b + 3)(a + 1) + ^(6 + 1)(6 + 2)(a + l)(a + 2) 
6 4 



1 

+6 



+ 



4 
1 

+ 2 



6 + l)(a + l)(a + 2)(a + 3) - ^(a + 1)(6 + 1)(6 + 2) - ^(a + l)(a + 2){b + 1) 

a + l)(a + 2)(a + 3) + ^(6 + 1)(6 + 2)(6 + 3) + ^(a + l)(a + 2)(6 + 1) 

6 2 

a + 1)(6 + 1)(6 + 2) + hb + l){b + 2){b + 3) - Ub + l)(a + l)(a + 2)(a + 3) 

6 6 

a + l)(a + 2)(6 + 1)(6 + 2) - ha + !)(&+ 1)(6 + 2)(6 + 3) + J(a + l)(a + 2)(6 + 1) 

6 2 

a + 1)(6 + 1)(6 + 2) + i(6 + 1)(6 + 2)(6 + 3) - ha + l)(a + 2)(6 + 1) 

6 2 

-^(a + l)(6+l)(6 + 2) 

J(a + l)(a + 2){a + 3) - ^(a + l)(a + 2)(6 + 1) + ^(a + 1)(6 + 1)(6 + 2) 



+i(6 + 1)(6 + 2)(6 + 3) = i(6 - a)(n + l)(n + 2), 



(167) 



where the sequence of terms in the ffist expression corresponds to the one in Eq. f ll63p . 
Therefore, expression fll63p can be rewritten as 



- E 



6 — a 



n 



:i68) 



a+b+2=n 

Adding to this expression the ffist term in Eq. ( 167|) . we obtain that the total derivative 
in Eq. ( j67ll can be rewritten in the following form: 



(169) 



a+b+2=n 



D Summation of subdiagrams with two external lines 



Here we describe, how subdiagrams with two external lines are split into groups, 
convenient for the calculation, and present results of this calculation. (Expressions for 
the left vertexes are omitted for simplicity.) 
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1. Subdiagrams with a chiral left end and an antichiral right end: 





After the substitution Y ^ O'^ the sum of these diagrams is written as 



9' 



(170) 



8(a2 + M2)3 " (92 + M2)2 16(92 + M2)' 

In particular, this means that in the massless case a sum of such diagrams is equal to 0. 
2. Subdiagrams with two chiral ends: 





These diagrams are given by 

0' ^'"^ 
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2M- 



2(a2 + M2)3 (a2 + M2)2 4(^2 + M2)' 



(171) 
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3. Subdiagrams with two antichiral ends: 





These diagrams are given by 



2M2 



2((92 + M2)3 V(52 + M2)2 4(^2 + M2) 2(^2 + M2)2 

4. Subdiagrams with an antichiral left end and a chiral right end: 



(172) 





These diagrams are given by 



4 2M^D^D^ 4 2M2 



(92 + M2)3 (92 + M2) 16(^2 + M2) ■ 



(173) 
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